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We show that the computational complexity of Riemann mappings can be bounded by the complex- 
ity needed to compute conformal mappings locally at boundary points. As a consequence we get 
first formally proven upper bounds for Schwarz-Christoffel mappings and, more generally, Riemann 
mappings of domains with piecewise analytic boundaries. 

1 Introduction 

The Riemann mapping theorem is probably one of the most fundamental theorems in complex analysis. 
Thus, not surprisingly, effectivization of this theorem was on the agenda ever since the theorem was 
proven. The first general constructive solution for Riemann mappings were given by R Koebe in 1910 
using the osculation method (see Ol). More constructivity results were added over the years based 
on different methods like potential theoretic methods, circle packings and random walks. Using the 
osculation method, R Hertling [ 4 ] gave the first exact characterisation of the effectiveness of the Riemann 
mapping based on type-2-theory of effectivity. Similar results were recently given for multiply connected 
domains (see [5]). 

Beside these effectiveness results it is known that the functor which maps simply connected domains 
to its Riemann mapping is jJP-hard (see e.g. [1]). Even for the Riemann mapping of a single domain, 
hardness results are known (see [9]). On the other hand most natural domains seems to have further prop- 
erties such as piecewise analyticity which seems to simplify the task of computing Riemann mappings of 
these domains. A prominent example of such a restricted class of domains are polygons and there exist 
seemingly fast algorithms to compute the Riemann mappings which are also called Schwarz-Christoffel 
mappings in these cases. A strict proof that theses algorithms are polynomial time bounded, however, 
is not known. Even worse, beside a few examples of domains for which the Riemann mapping is given 
explicitely, there is only one more general class known to have polynomial time computable Riemann 
mappings: Using Bergmann kernels one can show that domains with analytic boundaries can indeed 
be computed in polynomial time (see |9|). But even for polygons or domains with piecewise analytic 
domains this method breaks down and leads to exponential time bounded algorithms. Thus for these 
seemingly simple domains, the most efficient algorithms known so far are based on the random walk 
technique (see e.g. [1]) and need exponential time. 

In this paper we will present a more general class of domains with polynomial time computable 
Riemann mappings. Our results are based on a local to global principle. That is, we show, that if locally 
polynomial time computable Riemann mappings or, more generally, solutions to Dirichlet problems 
exist then also the Riemann mapping of the domain is polynomial time computable. As an immediate 
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consequence of this result we get polynomial time algorithms for polygons and, more generally, domains 
with piecewise analytic boundaries. 

To prove our main result we will use the well known connection between Riemann mappings and har- 
monic functions (see Section 2 below). Though there are well established solution methods for boundary 
value problems which work well in practice, it seems to be very hard to prove that these methods work 
indeed in polynomial time. The most promising method is probably the finite element method. Usually 
this method is using a more or less regular division of the domain, so that each of the subdomains (grid or 
triangulation) has roughly the size of the precision one wants to achieve. This clearly cannot give a poly- 
nomial time bounded algorithm. But even if we adopt this method and adjust the size of the subdomains 
similarly to the "layouts" which we will introduce in Section 3, there still is a problem which seems hai^d 
to overcome: To bound the time needed to solve the induced linear system one has to get some bound 
on the conditioning, i.e. absolute values of the determinants of the underlying linear systems. Although 
one knows that in the end everything will work out nicely there is no general method known to calculate 
this "end". Whereas such methods seems to be reachable for e.g. polygons, it is rather unlikely that we 
can compute bounds if we consider approximations to the boundary which can vary with the precision 
one wants to achieve. We will use a technique similar to the finite element method known in scientific 
computation, but we will neither compute a Unear system explicitely nor compute any bounds on the 
conditioning. Instead we will compute harmonic functions on suitable subdivisions and update these 
functions in several iterations. To analyse the time complexity of our method we will use a random walk 
technique rather than bounding determinants. 

We will end this section with a very short introduction to type-2-theory of effectivity. Afterwards in 
Section 2 we will prepare our main result which we will present in Section 3. Finally, in Section 4 we 
will discuss some consequences of the main theorem. 

Though we assume that the reader is familiar with the basic notions of type-2-theory as can be found 
e.g. in the textbook ifTTl . we have to adopt the model slightly to get a natural notion of complexity. To 
this end let £ be a finite alphabet and £*, £** denote the free monoid over £ and the class of functions 
£* — > £*, respectively. We say that a function / :C £** £** can be computed in time t iff there exists 
an oracle Turing machine M which computes f^{u) = v in at most t{\u\) steps whenever f{o){u) = v, 
i.e. the input of / is used as an oracle. Combining this with classical complexity and standard tupling 
functions we can extend this definition to any products of £* and £**. To define complexity on more 
interesting classes we will use representations and realizations: here a representation is a surjective 
functions v :C £* — > M or p :C £** — > M, where M is the class we want to represent. We will identify 
usual representations / :C N ^ M and g •.Q'L'" ^ M with representations v, p so that v(m) = f{\u\) 
and p(a) = g(a(£)a(0)a(00)...), respectively, where we assume G £. A realization F of a function 
f :Q Ml — M2 between two represented classes Mi, M2 with representations 5i and ^2, respectively, is 
a function F : dom(5i) — > dom(52) so that ^2 oF = / o 5i on the domain dom(/o 5i) of /o 5i. We say 
that / is computable in time t if there exists a realizer of F which is computable in time t. 

To simplify things further, we will next introduce standard representations and use these represen- 
tations later on without further mentioning. That is, we call a function computable in time t iff it is 
computable in time t with respect to the standard representations. 

We use the usual standard representations of the natural and dyadic numbers. Furthermore we use 
the signed digit representation for the reals. By identifying the complex numbers with pairs of reals we 
get immediately standard representations for N[/], Y[i] and C, respectively, where Y denotes the class of 
dyadics. Similarly, by identifying open balls ©^(z) = {z' € C | |z — z'| < r} with the tupel (r,z) we get a 
standard representation of the class B of such balls. By combinding this with a standard representation of 
the class of finite subsets of N we get our standard representation of the class B* of finite tuples of such 
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balls. Finally, we will need a suitable representation of harmonic functions. As we consider only a very 
restricted class of functions, i.e. all functions map some disc to [0; 1], the following representation omits 
any information on the modulus of continuity, which can be computed easily. For p G £** and m € £* 
let p||„ denote the function p||„ : £* ^ £* with p||i,(v) = p{uv). Then a name co € L** is mapped to a 
harmonic function / under our standard representation iff ft)||„ is a name of f{z) whenever m is a name of 
the dyadic z. Similar to balls we get a standard notation M* of all finite sequences of harmonic functions 
and finite sequences of other represented objects. Beside we will need piecewise harmonic functions. To 
simplify things we will in the sequel define for every s = {{fo,Bo),{fi,Bi), {fm^^m)) G (H x B)* a 
unique function h, : [j^Bi ^ [0; 1] by h,{z) := fi{z) iff z G BiXDj^^Bj. 

Beside the time needed to compute a certain function we will have to be careful about another re- 
source: The lookahead of a machine on some input a is the function /m(<7) : N — )• N which gives for 
given « G N the maximum length of the question asked by the machine to a on inputs of length n. We 
say that M has linear lookahead iff there exists a konstant k so that («) <k-n for all n and a. This 
resource is especially important for iterations of functions. Assume that we want to iterate a polynomial 
number of times a polynomial time computable function. If we are not carefuU about the lookahead this 
can easily lead to functions which are no longer polynomial time computable. If on the other hand we 
can guarantee that all these functions have Unear lookahead with a common constant then the function 
given by this iteration process is again polynomial time computable. 

2 Harmonic Functions, Dirichlet Problems and Riemann Mappings 

To simplify things we will in the sequel consider simply connected domains G which are bounded by a 
Jordan curve J so that D3/5 ^ G C D4/5. Any simply connected domain fulfilling the condition given 
in [4J can be transformed to such a domain by computable conformal mappings and, if the domain is 
bounded, even by dyadic polynomials. Let the class of these restricted domains be denoted by G. 

For G G G the classical Riemann mapping theorem states that for any z G G there exists a unique 
conformal mapping /(z) : G — > 'Bi{0) so that /(z) = and /'(z) > 0. Let us denote this map by for 
the moment. Obviously, any z so that can be computed in polynomial time on some neighborhood 
of is itself polynomial time computable. Furthermore, by using standard methods one can show that 
for z and z' in G, which are given via the standard representation of C, we can compute o {f^)^^ in 
polynomial time. Thus we can restrict ourselves to z = and we will write f*^ instead of for short. 

It is not hard to give examples of domains G G G so that cannot be computed in polynomial time 
throughout G because the modulus of continuity of this function tends to infinity for these examples 
when we approach the boundary of the domain. Thus, a polynomial time algorithm can only exist iff 
we compute only at values of a suitable distance to the boundary of G. Analyzing the mentioned 
examples shows that possible upper bounds on the modulus of continuity are given by its non-euclidean 
distance to z. To simplify things we will translate this to the euclidean distance d{-,-) and introduce the 
following notion of efficient computability inside G: 

Definition 1 We say that a function f is polynomial time computable inside G iff there exists a machine 
M and a polynomial p so that for a// « G N and z G so that d{z,dG) > 2^", M computs f{z) up to 
precision 2^" in time p{n). 

Notice that polynomial time computability on G does always imply polynomial time computability 
inside G for functions defined on G. Obviously, the converse does not hold. 

Using the above notion we can formulate the well established relation between polynomial time 
computable Riemann mappings and the solutions of Dirichlet problems (see iU) as follows: 
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Lemma 1 Let G G G and F : G ^ G be the function defined by F{z,w) = hy{w) for all z-,w ^G, where 
h^ denotes the harmonic function h, : G\Bg-2n(z) [0; 1] so that h^{v) = Ofor allv&dG and h^(y) = 1 
for all V G 5Dg-2n(z). If F is polynomial time computable inside G then the Riemann mapping f^ is 
polynomial time computable inside G. 

Now let G C D be a connected domain with regular boundary and (p : — [0; 1] be a continuous 
function on the boundary of G. Notice that indeed the domains we need in Lemma [T] have regular 
boundaries. Thus there exists a unique solution to the Dirichlet problem with boundary values given 
by (p. We will denote this solution by g^. It is well known that the Dirichlet problem can be solved 
easily on the unit disk D := Bi (0) via the Poisson kernel function P{w,z) = 1 /2n ■ (1 — |zP)/(|h' — zp). 
This function is harmonic on C \ {z.} for fixed z and furthermore we have that Jgj^(p{w)P{w,z)dw is a 
solution to the Dirichlet problem with boundary values given by (p. We will see in the next section that 
this solution can be efficiently computed if cp fulfills certain regularity conditions. The above integral 
can be generalized to arbitrary Jordan domains G by using the harmonic measure (Og{z), i.e. we have 
that g'^iz) = Jgc (p{x)d(OG{z){x) for all z€G. 

In our algorithm we will immitate the classical construction of solutions for the Dirichlet problem via 
subharmonic functions: Using the regularity of the boundary we can find for each boundary point a local 
barrier and thus a neighborhood of the boundary as well as a harmonic function on the common points 
of G and this neighborhood which coincede with (p on the boundary and is on all boundary points (of 
the neighborhood) which do not belong to the boundary of G. Fixing all other points in G to be one 
repeatedly modifies this function by taking some disk inside G and replacing the value of the function 
inside this disk by the values given by the integral Jgjjfh{w)P{w,z)dw where h denotes the values of the 
unmodified function on the boundary of the disk. Iterating this modification with different disks will in 
the end lead to the solution of the Dirichlet problem. We will need the following upper bound motivated 
by the above idea. 

Lemma 2 Let G (^ID) be a domain with regular boundary and ^ : dG ^ [0;l] be continuous. Then we 
have for all Jordan domains J QG and all intergrable functions *F : G — )• [0; 1] with *F < g'^: 

j ^ix)d(Ojiz)ix)<g'^{z) 

JdJ 

for all z(z J- 

The above idea motivated also the use of random walks to solve Dirichlet problems. This is an 
adopted version of Kakutanis theorem to our situation: 

Theorem 1 Let G(^Obe a domain, (pbe a continuous function on dG and let Rbe a continuous function 
/? : G — )• M w/f/i l/3-d{z,dG) < R{z) < 2/3 ■d{z,dG) for all z ^ G. Furthermore let the random process 
Xt (z) be defined by 

xoiz)=z, xtiz) =x^_i(z)+e2^'^' 

where Gq, 0i, ... is a sequence of independent random variables, uniformly distributed in [0; 1]. 
Then with probability 1, limt^ooXt{z) exists, limt^ocXt{z) G dG and E{(p{limt^ooXt{z)) = g^{z)- 

Discretizing this idea leads to the random walk algorithm in [IJ. To get the result to precision 2^" 
one has in each x, to take an exponential number of 0,. Thus clearly this cannot give polynomial time 
algorithms. On the other hand, however, the number t of iterations can be bounded by a polynomial to 
get precision 2^". We give this result by a modification of the above theorem. 
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where z' is some boundary point near x,_i [z), then 



\E{(p{xTiz)) 




2n 



where T = k2- n^ for a suitable constant k2 and (p{y) '■= Ofor all y ^ dG. 

The previous lemma will provide a lower bound by using the following lemma: 

Lemma 4 Let G, (p, R and Xt be given as in Theorem\J\ Furthermore let t and £ > be given. Then for 
any function *P : G — )• [0; 1] with *F(z) > E{(p{xt{z))) — £for all z £ G we have 



To get a starting point we thus have to find some neighborhood of the boundary of G and a (sub- 
)harmonic function on ^Fq on this neighborhood so that ^o{z) > E{(p{xi{z))) for all z- For the special 
domains which we will need by Lemma[T]this is nearly trivial: Let therefore for G S G, « G N and some 
Zo ^ G with d{zQ,dG) > 2^" by G"^ denote the domain G"^ = G\B^-2«(zo)- Furthermore let a covering 
Ni of the boundary dGbe given so that locally the Dirichlet problem can be computed efficiently on each 
of these neighborhoods. We will define this notion formally in the next section. Then we can choose *Po 
near the boundary, i.e. on the Ni, by = and furthermore choose on N"^ : — G'^^ nID'i/2-2 "(zo) to 

N" 

be the solution g^'° of the Dirichlet problem with boundary values given by (p where (p(5Dg-2/i(zo)) = 1 
and (p(5Di/2-2-"(zo)) = 0. Outside these sets we let *Po be *Po = 0. Then we compute the sequence ^Fq, 
^i*!, ^k^.ni where ^t+i is computed from *Pf by 



if z G A'^ for some A'^ — Ni ot N — ^"o' ^^^^ ^ denotes again a function with 1/3 • d(^z, dG) < ^(z) ^ 
2/3 • d{z, dG) for all z G G and ^, is given by ^P, = on dG, »Pf = 1 on dBi/2.2- (zo) and = on 
all other points. 

By Lemma 2, 3 and 4 ^i^^.ni is then an acceptable approximation of the corresponding Dirichlet 
problem which we will need to apply Lemma 1 . 




for all z& G. 




ifz^A^^'^UA^r-and 



3 The Main Theorem 



In this section we will show that the ideas developed in the last section lead indeed to a polynomial 
time algorithm for the Riemann mapping. Probably the easiest task of this proof is the characterization 
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of those G G G which allow to compute the function 7? of Theorem [T] effectively. To this end we have 
to compute the distance function to the boundary up to constant factors. Here the factor is somewhat 
arbitrary. We will call this class of domains polynomial time computable. Despite the fact that we have 
motivated the following definition by the ideas given in the previous section, the definition itself is quite 
natural and equivalent definitions have been used before to characterize computably efficient domains 
(see e.g. ifTOl). 

Definition 2 We call G G G polynomial time computable iff there exists a polynomial time computable 
function F : G ^ [0; 1] so that for all G we have (1 - 1/4) •J(z,5G) < F{z) < (1 + \/4) ■ d{z,dG). 

Let for the time beeing G G G be polynomial time computable and F its approximation function of 
the previous definition. Then the function /? : G — )• [0;1] with R{z) '■= 1/2 ■F{z) for all z G G is also 
polynomial time computable and we have indeed 

1/3 ■d{z,dG) < R{z) <2/3 ■d{z,dG) 

for all z ^ G. We even have more: Let r{z) '■= 1/4 • F(z). Then for every x G Or{z) (2) we have 

\/3-d{x,dG) <R{z) < 2/3 ■dix,dG). 

We will next discuss the computation of the integral /^jj^^ ^^^^^Wt{w)P{w,z)dw. It is well known that 
integration leads easily to [j-P hard problems and thus there is no polynomial time upper bound known. 
However, if we restrict ourselves to harmonic or analytic functions integration can be done in polynomial 
time (see e.g. [6]). The following result can be easily proven by such techniques. Let for a given ball 
B = Br{z) the ball D2r(z) be denoted by 2B. 

Lemma 5 Let G G G, d > be given. Furthermore assume that there exists a polynomial time com- 
putable function 5 : N — > (H X B)* so that with S{n) = {{fo,Bo), [fn^B,/)) we have 2B,- C G for all i 
and all functions f : B,- — )• [0; 1] can be extended to harmonic functions on 2Bifor all i. 

Then there exists a polynomially time bounded machine M with linear lookahead which computes for 
given n z € G and r €Y the integral 

/ hsu)iw)Piw,z)dw 

whenever Dr(z) is covered by the Bj and r > 2^'^'". 

Finally, we have to formally define what we mean by "locally efficient solutions of Dirichlet prob- 
lems". Roughly speaking, this means that we can cover the boundary by neighborhoods so that for each 
of these neighborhoods and any piecewise harmonic function on the boundaries of these neighborhoods 
one can compute the solution of the corresponding Dirichlet problem in polynomial time. 

Definition 3 An admissible neighborhood in a domain GQHis a tuple (z, r,D) where z G dG, r G Y and 
D(^Gis a simply connected domain so that Dsr^z) fl G C D and Osriz) CidGisa simple, connected curve. 
An admissible n-covering of such an admissible neighborhood is a finite sequence (Bq, ...,B,„^i) G B* 
so that 2Bj C Gfor i = l,...,m, the diameter of Bq and Bm+i <^f^ cit most 2^^'" and U/^i 5 dD\dG. 
Furthermore the core of an admissible neighborhood is the set Dr(z)- 

Let (jO : G ^ [0; 1] be given. We say that G C D admits locally polynomial time solutions of the 
Dirichlet problem with boundary values (p iff there exist polynomial time computable functions 5 : N — )• 
(D X Y X B*)* and U : dN^ xM* xC^C so that for all n there exists a sequence {zo,ro,Do), (zi,ri,Di), 

{ZrmfmjDm) of admissible neighborhoods in G so that 
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1. S{n) = {{zo,ro,Co),...,{Zm,rnt,C,„)), 

2. the cores of these admissible neighborhoods cover dG, 

3. the Q are admissible n-coverings of the {zi,ri,Di) and 

4. for all i gN and f = (/o, gW we have that g^' = U{n,i,f,-) wheneverCi = (Sq, 
q> = (pondGandp = h^(fy^Bo),...,{L+uBs+i)) ^" ^A' \ dG. 

S and U are then called solutions of the problem. The balls in Ci are called boundary nodes of the 
solution. 

With this preparations we can nearly implement the idea at the end of the previous section literally. 
Of course we can neither keep the values of (z) for all z nor can we update (z) to (z) for each 
z G G without any assumptions on We will construct % to be piecewise harmonic. Thus we can 
indeed keep track of any of the values of easily. To this end, however, we need some partition of G 
so that all *Pf are harmonic on each part of this partition. Furthermore there should not too many of these 
parts. We will call such partitions "layouts". The following definition will be used inside the proof of 
our main theorem. 

Definition 4 Let n gN be given, <p : G — > [0; 1] and U, S be locally polynomial time solutions of the 
Dirichlet problem with boundary values (p. Furthermore let S{n) = {{zo,ro,Co), {Zm-,rm-,Cm))- 

We say that (So, ■■■,Bs) G B* is an n-layout for this solution iffG\ U,ID'3r; C \JjBj and for each Bj 
we have that 2Bj C G. .v is called the size of the n-layout. 

Now let G be polynomial time computable. We say that G admits polynomial time computable layouts 
iff there exists a polynomial time computable function V :N ^M* so that V{n) is an n-layout and for 
each Bir{z) in V{n) we have r = r{z) = 1 /4F{z). 

The last step we need to prove our main theorem is a proof that polynomial time computable layouts 
exists. We will give the main idea of such a proof. A formal proof will be given in the full version of this 
paper. 

Lemma 6 Let G C.ID be polynomial time computable. Furthermore /e? <p : G — > [0; 1] be given so that 
G admits locally polynomial time solutions U and S of the Dirichlet problem with boundary values (p. 
Then G admits also polynomial time layouts for this solution. 

Proof (Idea): 

First notice, that the construction of layouts is simple. We can start with a point z inside G, compute 
the correct radius of the corresponding disc, i.e. 1/4 ,F(z), and continue with some point on the boundary 
of this disc, while checking wether G is already covered by the neighborhoods given by U and the 
constructed discs. It is thus not hard to prove the following claim: 

Claim 1 G admits polynomial time layouts for the solution iff there exists a polynomial p so that for 
each n there exists an n-layout of size p{n) for the solution. 

To show that such n-solutions of polynomial size always exist we proceed as follows: Let be 
some connected component of G\|J;3B;. where the B,'s are the cores of the corresponding admissible 
neighborhoods. W.l.o.g. we can assume that there exists some d > (common for all n) so that all 
boundary nodes and all B,'s have diameter at least 2~'^". Now let z G dK be given so that d{z, Ui^O is 
minimal and add the disc ^r{z) (^) ^ layout. Furthermore mark the core S; by level 0. Next proceed 
with Kq replaced by Ki = K \I])2r[z){z) and U/^i replaced by Dr(z)(^) UU;^;- Notice, that any in this 
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way added disc has marked some core. Generally, if j is minimal so that z belongs to dKj we mark the 
corresponding core by level / 

It is not hard to show that each core can be marked at most 671 times at each level. Furthermore the 
maximum level by which a core can be marked is bounded by logj § 2'^ " because if we go from level k 
to ^ + 1 the diameters of the considered discs grow at least by a factor (1 + ^3/4) as can be seen by 
elementary geometry. 

Thus in the end we get a layout of size c m-n for a. suitable constant c (which does not depend on 
n) where m denotes the number of cores. As this number is polynomially bounded by the definition of 
locally polynomial time solutions of Dirichlet problems, we get altogether a polynomial upper bound on 
the number of discs in the layout. 

□ 



Now our main result follows easily by the idea given at the end of the previous section: 

Theorem 2 Let G G G be a polynomial time computable simply connected domain which admits locally 
polynomial time solutions of the Dirichlet problem with boundary values = 0. Then, for any z G, the 
corresponding Riemann mapping ff is polynomial time computable inside G. 

Proof (Idea): 

By Lemma [T] the theorem is proven once we have shown that the Dirichlet problem on G"^ with 
boundary values = and = 1 can be proven in time p{n) where is a suitable polynomial. 

First notice that G admits locally polynomial time solutions of the Dirichlet problem iff G"^ admits 
such solutions because the boundary componend (zo) admits clearly a constant number of admis- 
sible neighborhoods and corresponding solutions. By Lemma[6]we can therefore compute a polynomial 
time layout for this solution. 

Let us call any of the boundary or inner nodes just nodes and let B, be the cores of the solution. We 
start with exactly the *Po given at the end of the previous section. Then we update every node D = B,-(z) 
as follows: If D n IJ,- 3B, ^ then we have by definition of admissible neighborhoods that 2D C 5Bj 
for some Bj and thus we compute the new values on D by computing the solution of the local Dirichlet 
problem on the corresponding neighborhood. Otherwise we compute the new values on D by integration 
(Lemma [5ll. 

Let us call the procedure of updating all nodes a round. Each round can be computed in polynomial 
time. Thus to proof our main theorem it remains to show that we only need a polynomial number of 
rounds to approximate the solution of the Dirichlet problem to precision 2^". However, at the end of 
each round we get some piecewise harmonic function on G and a simple inductive argument shows that 
the piecewise harmonic function f we get after t rounds is larger than the function *Ff which we have 
introduced at the end of the previous section. Especially we get ft{z) > E{(p{xt{z))) where (p is the 
corresponding boundary value (=0 ondG and = 1 on the other boundary component). Thus we get that 
after T = k2- n^ rounds we have that fr is larger than the solution minus 2^^". By Lemma[2]we have on 
the other hand that for any t the function /, must be smaller than the solution. Thus after T rounds we 
get a correct approximation of the solution. 

As the number of nodes is polynomially bounded and any update can be done in polynomial time, 
we iterate polynomial times a polynomial time algorithm with linear lookahead. That is the overall 
algorithm is indeed polynomial time bounded. 

□ 
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The following well known result shows that it suffices to consider approximations of the boundary of 
G in the above construction. Thus the above result although holds iff locally polynomial time solutions 
of the Dirichlet problem exists for approximations of G, where these approximations may depend on the 
precision one wants to achieve. We omit the details. 

Lemma 7 Let ji and 72 be two closed Jordan curves with d{yi , 72) < 2^^". Furthermore let hi and h2 
be two continuous functions on 71 and ji, respectively, with the following continuity property: ifx\ G 71 
and X2 € 72 and \x\ — X2I < 2^^" then \ui{xi) —U2{x2)\ < 2^". Then the solutions Uf,^ and ui^^ of the 
corresponding Dirichlet problem fulfills 

\uh,{z)-Uh^iz)\<2-"+^ 
for all z inside 71 and 72 with d{z, 7i ) > 2^" and d{z, 72) > 2^". 

4 Discussion 

We end this paper with a short list of applications of and some remarks on our main theorem. 

With our main theorem it is now easy to prove that Schwart-Christoffel mappings are polynomial 
time computable where we do not even have to mention the parameter problem. To define admissible 
neighborhoods we have simply to use the fact that mappings for triangles can be computed in polynomial 
time because here the parameter problem is trivial (see e.g. |3|). Thus we can choose simply triangles 
as local neighborhoods and use the Schwarz-Christoffel mappings and the relection principle to solve 
Dirichlet problems locally: 

Corollary 1 Let G be a polygon so that any vertex is polynomial time computable. Then for any z€ G 
the Riemann mapping f^ is polynomial time computable. 

A similar idea shows that we can extend this result to piecewise analytic boundaries. Here we use the 
fact that we can map the analytic boundary locally one to one onto some intervall. To get an admissible 
neighborhood we therefore just have to use this map, take a small triangle and map it back by the inverse 
of the first mapping. Thus by applying again Schwarz-Christoffel mappings for triangles we can solve 
locally Dirichlet problems. 

Corollary 2 Let G be a simply connected domain with piecewise analytic, polynomial time computable 
boundary. Then for any z £G the Riemann mapping f^ is polynomial time computable. 

Notice that this second result cannot be extended to domains where the boundary is approximated 
by piecewise analytic boundaries which may differ for different precisions. To this end, more results on 
what kinds of approximation could lead to polynomial time algorithms would be interesting. 

On the other hand such results can be easily shown for polygons where the number of vertices is 
bounded polynomially by n and the vertices can be computed uniformly in polynomial time. 

Combining the previous result with the result |9| we get an improved lower bound on the complexity 
of Riemann mappings for single domains. Details can be found in the original paper. 

Corollary 3 /f jjf ^sei P then there exists a polynomial time computable simply connected domain G so 
that for any computable z ^G we have that f^ is not even polynomial time computable on a compact 
subset ofG. 
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